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Is the common envelope ejection efficiency a function of 
the binary parameters? 
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ABSTRACT 

We reconstruct the common envelope (CE) phase for the current sample of observed 
white dwarf-main sequence post-common envelope binaries (PCEBs). We apply multi- 
regression analysis in order to investigate whether correlations exist between the CE 
ejection efficiencies, cvce inferred from the sample, and the binary parameters: white 
dwarf mass, secondary mass, orbital period at the point the CE commences, or the 
orbital period immediately after the CE phase. We do this with and without consid- 
eration for the internal energy of the progenitor primary giants' envelope. Our fits 
should pave the first steps towards an observationally motivated recipe for calculating 
acE using the binary parameters at the start of the CE phase, which will be useful 
for population synthesis calculations or models of compact binary evolution. If we do 
consider the internal energy of the giants' envelope, we find a statistically significant 
correlation between acE and the white dwarf mass. If we do not, a correlation is found 
between cvce and the orbital period at the point the CE phase commences. Further- 
more, if the internal energy of the progenitor primary envelope is taken into account, 
then the CE ejection efficiencies are within the canonical range < «ce < 1, although 
PCEBs with brown dwarf secondaries still require «ce ^ 1- 
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1 INTRODUCTION 

The common envelope (CE) phase is a key formation process 
of all compact binary systems, such as cataclysmic variables 
(CVs) and double white dwarf binaries. Such systems typ- 
ically have orbital separations of a few solar radii, and yet 
require orbital separations of between approximately 10 and 
1000 to accommodate the giant progenitor primary star. 

IPaczvnskH (|l976t ) first suggested that the CE phase is 
responsible for removing large amounts of orbital energy and 
angular momentum from the progenitor system, causing a 
significant reduction in the orbital sep aration between the 
two stellar com ponents (for reviews see llben fc Livioi ri993: 
IWebbink|[20ul ). 

In the pre-CE phase of evolution, the initially more mas- 
sive stellar component (which we henceforth denote as the 
primary) evolves off the main sequence first. Depending on 
the orbital separation of the binary, the primary will fill 
its Roche lobe on either the red giant or asymptotic giant 
branch (AGB) and initiate mass transfer. If the giant pri- 
mary possesses a deep convective envelope [i.e. the convec- 
tive envelope has a mass of more than approximately 5 per 
cent of the giant's mass; iHiellming fc Webbinkl §987)], the 



giant will expand in response to rapid mass loss. As a result, 
the giant's radius expands relative to its Roche lobe radius, 
increasing the mass transfer rate. As a consequence of this 
run-away situation, mass transfer commences on a dynam- 
ical timescale. The companion main-sequence star (hence- 
forth the secondary) cannot incorporate this material into 
its structure quickly enough and therefore expands to fill its 
own Roche lobe. The envelope eventually engulfs both the 
core of the primary and the main-sequence secondary. 

During the spiral-in phase, if enough orbital energy is 
imparted to the CE before the stellar components merge, 
then the CE may be ejected from the system leaving the core 
of the primary (now the white dwarf) and the secondary star 
at a greatly reduced separation. We term such systems post- 
common envelope binaries (PCEBs). In the present study, 
we consider PCEBs which have a white dwarf primary com- 
ponent, and a main sequence secondary companion. 

Modelling the CE phase presents a major compu- 
tational challenge, and can only be adequately accom- 
plished vi a three-dimensional (3 D) hydrodynamical simula- 
tions (e.g. ISandauist et al.ll2000l ). Even this approach, how- 
ever, still cannot cope with the large dynamic range of time 
and length scales involved during the CE phase. 



2 P. J. Davis et al. 



Clearly, such hydrodynamical simulations are too com- 
putationally intensive to be included in full binary or pop- 
ulation synthesis codes. Instead, such codes resort to pa- 
rameterisations of the CE phase. One such parametrization 
describes the CE phase in terms of a simple energy budget 
argument. A fraction qce of the orbital energy released as 
the binary orbit tightens, A_E or b, is available to unbind the 
giant's envelope from the core. Hence, if the chan ge in the 



envelope's binding energy is A_Ebmd, we have (e.g. de Kooll 
Il992l ; Ide Kool fc Ritterll 19931 ; IWillems fc Kolbll2004l ) 



AE bin 



(1) 



The efficiency oce is a free parameter with values < 
Qce < 1 (although values of qce > 1 are discussed). How- 
ever, the value of qce is poorly constrained as a result of 
our equally poor understanding of the physics underlying 
the CE phase. 

An alternative formulation in terms of the angu- 
lar momentum budg et of the binary was suggested by 
iNelemans et al.l (|2000l ) from their investigation into the for- 
mation of double white dwarf binaries, which they assumed 
occurred from two CE phases. They suggested that the rela- 
tive change in the binary's total angular momentum, A J/ J, 
and the relative change in the binary's total mass, AM/M, 
during the CE phase are related by 



AJ AM 

= 7CE" 



J 



M 



(2) 



where 7ce is the specific angular momentum removed from 
the binary by the ejected CE in units of the binar y 's ini - 
tial specific angular momentum. INelemans fc Tout] (|2005l ) 
found that a value of 1.5 < 7ce < 1-75 could account for all 
observed PCEBs. 

This alternative description was prompted by the fact 
that for double white dwarf binaries a value of qce < 
was needed to describe the first CE phase, which is clearly 
unphysical. This indicates that the orbital energy increases 
during the first phase of mass transfer, with a corr esponding 
increa se in the orbital separation. In light of this, IWebbinkl 
(2008) suggested that the first phase of mass transfer is 
quasi-conservative, and does not give rise to a CE phase. 

However, the angular momentum budget approach pre- 
dicts PCEB local space densities that are appr oximately a 
facto r of ten larger than observed estimates (|Davis et al.l 
l20ld ). It also predicts that the number of systems increases 
towards longer orbital periods, with a maximum number 
of systems at app roximately 1000 d jMaxted et al.l l2007al ; 
iDavis et al.ll20101 ). This is contrary to observations, which 
show that there is a maxi mum number of PCEBs with or - 
bital periods of about 1 d (|Rebassa-Mansergas et al"] [2008). 
Furthermore, the angular momentum budget approach has 
far less predictive power than the energy budget description; 
a tightly constrained value of 7ce cannot constrain the pa- 
rameters of the possible progenitors of observed systems, 
while th e energy budget approach is more prom ising in this 
respect (|Beer et al.ll2007l ; IZorotovic et al.ll20icl ). 

The energy budget approach therefore appears to be the 
more favourable approach. However, a few major unresolved 
issues have yet to be overcome if we are to make significant 
progress in understanding and modelling the CE phase. 

The first issue is that we require qce > 1 to account 
for some observed systems such as the double white dwarf 



binary PC 1115+166 l|Maxted et all I2002T) and the white 
dwarf-main sequence PCEB IK Peg jDavis et~ai1l2010r ). This 
indicates that an energy source in addition to gravitational 
energy is being exploited during the CE ejection process 
that is not accounted for in the standard energy budget for- 
mulation such as thermal energy and recombinat ion energy 
of ionized material within the giant's env elope (|Han et al.l 
2002; Dewi fc Taurisllioool ; IWeb bink 20o3). 

The second issue is that previous studies into the 
formation and evoluti o n of white dwarf bin aries (e.g. 
Ide Kool fc Ritterl 1 19931 ; I Willems fc Kolbl I2004T I have as- 
sumed that qce is a global constant, i.e. qce is the same for 
all progenitor systems, irrespective of their parameters upon 
entering the CE phase. This may be a so mewhat naive as- 
sumption. Indeed. pTerman fc TaamI (1 19961 ) suggest that qce 
may depend upon the internal structure of the Roche lobe 
filling primary progenitor. 

In this spirit, iPolitano fc Weilerl <|2007h calculated the 
present-day population of PCEBs and CVs if qce is a func- 
tion of the secondary mass. They considered both a power- 
law dependence (i.e. qce oc M2 n , where n is some power), 
and a dependence of qce on a cut-off mass. If the secondary 
mass was below this cut-off limit, then a merger between 
the two stellar components is assumed to be unavoidable 
(formally qce = in this case). 

As lDavis et al l (|2010f > found, however, these formalisms 
tend to underestimate qce for PCEBs with low-mass sec- 
ondaries (<0.2 Mq). For example even though a n — 2 
dependence could account for IK Peg, it predicts that sys- 
tems with M2 < 0.4 Mq cannot survive the CE phase. This 
is in conflict with observations that show that PCEBs with 
secondary masses as low as M2 ~ 0.05 Mq do exist. 

During the preparation of this paper, we became aware 
of sim ilar studies being c a rried out by Zorotovic et al.l 
j2010T ) and iDe Marco etail i|201ll ). IZorotovic et all (|2010h 
applied the CE recon struction method described by 
INelemans fc Tout! 120051 ) to the observed sample of PCEBs. 
They investigated whether qce depended on either the sec- 
ondary mass or the orbital period immediately after the CE 
phase. They concluded that there was no dependence in ei- 
ther case. However, this conclusion was based on a qualita- 
tive, 'by-eye', analysis. 

157 Marco et alj (|201ll ). on the other hand, found evi- 
dence that qce increases with increasing mass of the pro- 
genitor primary star, but decreases with increasing mass of 
the secondary. However, they only considered PCEBs which 
underwent negligible angular momentum loss since emerg- 
ing from the CE phase, which limited their sample size to 
30 systems. 

Our investigation is similar to these studies, except that 
we investigate whether qce can be treated as a first-order 
function of a combination of the binary parameters - white 
dwarf mass, progenitor primary mass and orbital period at 
the start of the CE phase, secondary mass and post-CE or- 
bital period - by performing multi-regression analysis. Our 
method for reconstructing the CE phase f or the observed 
PCEB sample is similar to that used by IZorotovic et alj 
|2010l) . but has been independently developed. The aim is 
to determine whether correlations exist between qce and 
the binary parameters. The resulting fits will then be useful 
recipes for binary evolution and population synthesis codes. 
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Such fits may also be useful to CE theoreticians in order to 
inform and constrain their models. 

Our PCEB samp le consists of systems contained in the 
iRitter fc Kolbl |2003t ) catalogue, Edition 7.14 (2010) and 
the 34 new PCEBs detected from the Sloan Digita l Sky 
Survey (SDSS), in contrast to |Pe Marco et al.l l|201lf) . For 
an in-depth review of this sample, we refer the reader to 
IZorotovic et~all l|2010l) . 

The structure of the paper is as follows. In Section 2 
we describe our method for reconstructing the CE phase 
for the observed sample of PCEBs. Our results are then 
given in Section 3, which are then discussed in Section 4. 
W e then compare our wo rk with similar studies carr ied out 
bv lZorotovic et al.l l^Old ) and lDe Marco et all (^OHl - ) in Sec- 
tion 5. Finally, our conclusions are given in Section 6. 



2 COMPUTATIONAL METHOD 

In the following subsections, we describe the energy budget 
equation describing the CE phase in more detail. We then 
explain how we can calculate the orbital separation of the 
system immediately after the CE phase from the cooling age 
of the white dwarf and the assumed angular momentum loss 
rate. Finally, we discuss how we use our population synthesis 
code BiSEPS, to evolve the possible progenitors of a given 
observed system through the CE phase. We then use the 
energy budget equation (c.f. eqn. [TJ to solve for qce- 



2.1 Calculating qce 

We calculate the binding energy of the primary's envelope, 
-Ebind, using: 

GMlMenv 



Ebin 



APi 



(3) 



where Mi is the mass of the progenitor primary, M env is the 
envelope mass and Pi is the radius of the primary, which is 
also equal to its Roche lobe radius at the start of the CE 
phase. The parameter A is the ratio between the approxi- 
mate expression of the gravitational binding energy given 
by eqn. @ and the exact expression: 



Ml 



GMi (r) 



dm, 



(4) 



where M c is the core mass (i.e. the mass of the proto-white 
dwarf) and Mi(r) is the mass enclosed within a radius r. 
Alternatively, if we consider the internal energy of the enve- 
lope, we have 

" M i r„\ r M i 



-Ebind 



GMi(r) 



dm + a th 



M c 



Udm, 



(5) 



where U is the energy per unit mass, and a t h is the frac- 
tion of thermal energy w hich is used to unbind the enve- 
lope ijDewi fc Taurisll2000h . We calculate qce for observed 
systems by adopting either eqn. ([4J (A = A g ) or eq. ([5]) 
(A = Ab) to calculate A in eqn. Q. Furthermore, we follow 
iDewi fc Taurisl (|2000h and assume that a t h = 1 for this case. 

For a binary system which has an initial orbital sepa- 
ration ylcE,i immediately before the CE phase (i.e. at the 
point when the primary giant just fills its Roche lobe), and 
an orbital separation Acej immediately after the CE phase 



(i.e. the point when the PCEB emerges from the CE phase), 
the change in orbital energy is given by: 

GM C M 2 GM1M2 



AP or b = 



(6) 



2^1-CE,f 2^4cE,i 

where M2 is the mass of the secondary star. 

By combining eqns. {TJ, ((3| and ((6| we can solve for 
qce- This gives, 



2M env Mi 



nM 2 A 1 M c A C e,i - Mi A CEA 



(7) 



where ri = Ri/AcE,i is the radius of the primary star in 
units of the initial orbital separation. 

2.2 Calculating Aces for observed systems 

Generally, the observed orbital separation of a PCEB is 
not the same as the orbital separation immediately after 
the CE phase. This is because the PCEB would have un- 
dergone orbital period evolution due to angular momen- 
tum losses via gravitational radiation or magnetic brak- 
ing since emerging from the CE phase. We therefore follow 
ISchreiber fc Gansickel l|2003h to calculate the orbital period 
of our observed sample of PCEBs immediately after their 
emergence from the CE phase, Pce, using the cooling time 
of the white dwarf, i coo i. 

We apply the di s rupted magnetic braki ng paradigm 
i|Spruit fc Ritterlll983l : iRappaport et al.lll983l ) and assume 
that magnetic braking is ineffective for secondaries that have 
masses less than the fully convective mass limit for mai n 
sequence stars, M con v,MS = 0.35 M l|Hurlev et al.ll2002l ). 
Thus, for observed PCEBs with M2 < A/ conVi MS, we assume 
that the systems are driven purely by gravitational wave ra- 
diatio n. We therefore calculate Pc e according to equation 
(8) of lSchreiber fc Gansickel (|2003l ): 



D 8/3 
CE 



256G 2/3 (27r) 8 / 3 t cool 
5c 5 

M1M2 8/3 

(Mi + M 2 y/ :i orb ' 



(8) 



The evolution of observed systems with M2 > M conv ,MS 
will be driven by a combination of gravitational radiation 
and magnetic braking. In this investi gation, we adopt th e 
magnetic braking formalism given by iHurlev et al. (2002). 
However, we can neglect the contribution due to gravita- 
tional radiation because the associated angular momentum 
loss rate is much less than that associated with magnetic 
braking. To calculate Pce we therefore apply 



,10/3 



10(2Tv) 10/3 rj h t cool R^Mc. 
1.72 x 10 15 G 2 / 3 
(Mi + M 2 ) 1/3 



+ P. 



10/3 



(9) 



where Pce, Porb and t coo \ are expressed in yr, and G is ex- 
pressed in Rq Mq 1 yr -2 . We also include a normalisation 
factor ?7h = 0.17 in order to obtain an angular momentum 
loss rate that is appropriate for the observed w idth and loca- 
tion of the CV period gap l|Davis et al.ll2008l ). We estimate 
the convective envelope mass of the secondary star as 



M c , 



M . 







M 2 > 1.25 M 







0.35 (i^ 



M 2 * 



0.35 < M 2 / M < 1.25 



(10) 
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l|Hurlev et aJj|20o3 ). Note that equation (JTOj) calculates the 
convective envelope mass for a star on the zero-age main 
sequence (ZAMS). In reality, the mass of the convective en- 
velope is a function of the fraction of the main sequence 
lifetim e that the star has passed through (e.g. iHurlev et al.l 
2000), thus making equation (|10[1 time dependent and the 
solution to Pce in equation ([9]) non-trivial. 




2.3 The modified BiSEPS code 

We use a modified ver sion of the BiSEPS code 
l|Willems fc Kolbll2002l . 120041 ) to evolve all possible progeni- 
tors of each observed PCEB through the CE phase, and to 
calculate «ce for each of them. 

The BiSE PS code emplo ys the single star evolution for- 
mulae bv lHurlev et alj 12000}) and a binar y evolu tion scheme 
based on that described bv lHurlev et al] |2002ft . We evolve 
a large number of binary systems initially consisting of two 
ZAMS stellar components. The stars are assumed to have a 
Population I chemical composition and the orbits are circu- 
lar at all times. The initial secondary mass is equal to the 
observed secondary mass of the observed PCEB we are con- 
sidering. The initial primary masses and orbital periods are 
obtained from a two-dimensional grid consisting of 1000 log- 
arithmically spaced points in both primary mass and orbital 
period. The initial primary masses are in the range of 0.1 
to 20 Mq, while the initial orbital periods are in the range 
of 0.1 to 100 000 d. Hence, we evolve 1 x 10 6 binaries for a 
maximum evolution time of 10 Gyr. For symmetry reasons 
only systems where M2 > M\ are evolved. 

Prior to the onset of the CE phase, the primary star will 
lose mass via stellar winds. For stars on and after the main 
sequence, ma ss loss rates are calcul a ted u sing the formalism 
described by iKudritzki fc Reimersl (|l978ft . The super- wind 
phase on t he AGB is calculated usin g the prescription de- 
scribed bv lVassiliadis fc Wood! (| 1993ft . 

If a binary configuration undergoes a CE phase, and the 
mass of the primary core is equal to the observed white dwarf 
mass within the observed uncertainty then we can use eqn. 
to calculate qce for that system. For each progenitor, 
we calculate A g or Ab by linear inte rpolating between values 
tabulated bv lDewi fc Taurij (|2000ft . We extended this table 
using full stellar models calculated by the Eggleton code 
(provided by Marc van der Sluys, private comm unication) 
and the EZ s tellar evolut i on co de (|Paxtonll2004ft . We refer 
the reader to iDavis et al.1 l|2010ft for further details. 

For a given observed PCEB, we now have a range of 
possible values of cvce, each corresponding to a possible pro- 
genitor system. A value of qce is given a weighting equal to 
the formation probability of the associated progenitor sys- 
tem. The formation probability for a progenitor can be found 
from the initial mass function (IMF), /(Mi,;), where Mi,i is 
the ZAMS primary mass (see eqn. IA1J1 . and from the initial 
orbital period distribution H(P rb,i), which can be found 
by using the initial orbital separation distribution (IOSD), 
h(ai) (see eqn. IA2[I . via Kepler's Third Law. The quanti- 




Figure 1. Distributions in the possible primary progenitor mass 
at the start of the CE phase (Mi, bottom panel), orbital pe- 
riod at the start of the CE phase, Pqe h (middle panel) and 
the CE ejection efficiency, oce (where we use A = A g ; top 
panel) for SDSS2216+0102 (red), SDSS1548+4057 (blue) and for 
SDSS0303-0054 (green). The vertical lines indicate the mean val- 
ues for each histogram. 

ties P or b,i and a; are the initial orbital period and orbital 
separation at the ZAMS stage respectively. 

If t he internal energy of th e envelope is considered (i.e. 
A = Ah). lDewi fc Taurisl (|2000f ) find that the binding energy 
of the envelope is positive if the radius of the primary gi- 
ant star is sufficiently large, formally indicated by a value 
of Ab < 0. This is a result of the fact that the magnitude of 
the radiation pressure exceeds the gravitational force. As a 
result, the prescription given by eqn. Q breaks down, and 
we cannot follow the evolution of such progenitors through 
the CE phase using the energy budget formalism (the impli- 
cations of Ab < is discussed in more detail in Section 3}. 
Therefore, progenitor binaries of a given PCEB with such 
primary giants are discounted in the aforementioned weight- 
ing calculation. 

For the secondary mass, we just consider a flat initial 
mass ratio distribution, g(qi) = 1 (see eqn. IA3fl . where qi = 
where M% t i is the mass of the ZAMS secondary. 
Hence, the probability that a PCEB formed from a CE phase 
with an ejection efficiency in the range Aqce is proportional 
to /(Mi,i)A(a I )ff(®). 

Fig. [T] shows distributions of possible primary progeni- 
tor mass at the start of the CE phase (Mi, bottom panel), 
orbital period at the onset of the CE phase (PcE,i, middle 
panel) and the CE ejection efficiency («ce, top panel, where 
we use A = A g ) for SDSS2216+0102 (red), SDSS1548+4057 



1 We can assume that the core of the primary giant is equal to the 
white dwarf mass because the increase in the core's mass during 
the CE phase will be negligible. 



2 Note that it does not matter how we determine the formation 
probability of the secondary star, as it will be the same for all 
progenitors of a given PCEB. 
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y = M„ D y = M, y = M 2 y = P CE y = P, 




Figure 2. Correlation matrix plot for the considered binary parameters. Rows correspond to independent variables, while columns 
correspond to dependent variables. Errors and numbers are omitted for clarity. 



(blue) and SDSS0303-0054 (green). The vertical lines indi- 
cate the mean value for each histogram. 

The white dwarf masses for SDSS2216+0102, 
SDSS 1548+4057 and SDSS0303-0054 are 0.400 ±0.060 M , 
0.646 ± 0.032 M and 0.912 ± 0.034 M Q respectively. We 
can see from the bottom and middle panels that the average 
primary mass and orbital period immediately before the 
CE phase increase for increasing mass of the white dwarf 
descendant. 

The orbital periods of these systems immediately after 
the CE phase are approximately the same (« 0.2 d). There- 
fore, the typical progenitor of SDSS0303-0054 would need to 
spiral-in further during the CE phase than SDSS2216+0102 
in order to reach the same post-CE orbital period. There- 
fore, for the former system, more orbital energy would be 



needed to eject the CE from the system than for the lat- 
ter system. Hence, the typical ejection efficiency will be less 
efficient for SDSS0303-0054 than for SDSS2216+0102. This 
is illustrated by the histograms in the top panel. The av- 
erage orbital separation and primary mass before the CE 
phase for SDSS1548+4057 lies between the average values 
for SDSS2216+0102 and SDSS0303-0054. Therefore, the av- 
erage ejection efficiency for this system also lies between 
SDSS2216+0102 and SDSS0303-0054. 



2.4 Statistical Approach 

Our aim in this investigation is to apply a chi-squared multi- 
regression analysis (via the Levenberg-Marquardt method) 
to obtain fits of qce for two cases; we do this by first deter- 
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Figure 3. Our reconstructed values of qqe versus the following binary parameters, where we take A = A g : white dwarf mass, M^oi 
primary mass at the moment the CE phase commences, Mi; secondary mass, M2; orbital period at the moment the CE ends, Pce; 
logio of the orbital period at the start of the CE phase, PcE,i- The solid green lines indicate linear fits to the data. 



Table 1. A summary of the fit parameters eo and ei for the data in each 
panel of Fig. [3] where we consider A = A g . Also shown for each panel are 
the values of T and the probability of exceeding T, ,^(P). Finally, an 
indication of whether we can reject the null hypothesis is included 

(see text for further details). 



Q 


£0 


ei 


T 


Z?(T) 


Reject ? a 


A/wd 


-0.13±0.05 


-0.48±0.27 


3.2 


1 x ict 1 


No 


Mi 


0.06±0.09 


-0.27±0.18 


2.2 


1 x 10- 1 


No 


M 2 


-0.26±0.08 


-0.42±0.15 


7.4 


9 x icr 3 


Yes 


Pce 


0.07±0.04 


0.39±0.07 


29.5 


< 10- 6 


Yes 


PcE,i 


0.94±0.33 


-0.34±0.11 


9.1 


4xl0" 3 


Yes 



a ^crit = 4.00 for a 0.05 significance level. 



mining whether correlations exist between acB and the bi- 
nary progenitor parameters at the point that the CE phase 
commences (i.e. the moment the primary progenitor star 
just fills its Roche lobe), which will be useful for popula- 
tion synthesis calculations. In the second case, we determine 
whether correlations exist between qce and the (observ- 
able) PCEB parameters, which may be a useful diagnostic 
to probe CE physics. 

We proceed as follows. We fit our reconstructed values 
of acE according to the equation 



log 10 Q C E = A) + ^2 Pilog 10 Qi, (11) 

i=l 

where /3o...m are constants. If we are considering 
the progenitor binary parameters only, then Q £ 
{-PcE.i, M2, Mwd, Mi}, with m = 4. If we consider the 
PCEB parameters, then Q £ {-Pce, M2, Mwd}, with m = 3. 

Strictly speaking, we should only apply the Levenberg- 
Marquardt chi-squared method if the errors are normally 
distributed. As we can see from Fig. [TJ this is not the case. 
We may therefore under-estimate the uncertainties in our 
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Table 2. A summary of our fit parameters, obtained from a multi-regression analysis of the values of oqe (where 
X = A g ) as a function of two or more binary parameters. Shown are the coefficients 00. ..m, the F-statistic, J 7 , and the 
probability of exceeding J 7 , 2P{T). We also indicate whether we can reject the null hypothesis J^q ■ 



Variables 


00 


01 


02 03 04 


T 


»(T) 


Reject J?* ? 








Progenitor parameters 








-PcE.i + -^WD 


2.35 


-0.78±0.25 


1.09±0.56 


3.6 


6xl0~ 2 


No 


+M 2 


1.94 


-0.67±0.27 


1.04±0.56 -0.19±0.19 


1.1 


3 x 10- 1 


No 


+Mi 


1.88 


-0.66±0.30 


0.97±1.38 -0.20±0.20 0.04±0.72 


0.003 


9.6 x 10" 1 


No 








Observable (PCEB) parameters 








PCE + Af WD 


-0.03 


0.20±0.03 


-0.85±0.21 


16.6 


2 x 10~ 4 


Yes 


+M 2 


-0.31 


0.28±0.02 


-0.95±0.13 -0.17±0.09 


38.4 


< 1 x lO" 6 


Yes 



10. a 7 




1 

a CE (mod.) 



Figure 4. Our reconstructed values of oce (y-axis, acE(obs.)) 
versus values of q?ce (x-axis, acE(mod-)) which have been calcu- 
lated from eqn. H14I I. Also shown is the value of the intercept and 
gradient of the linear fit through the data, which is shown as the 
solid black line. 




1 

a CE (mod.) 



Figure 5. Same as Fig. [4] except acE(mod.) has been calculated 
from eqn. 1151 . 



fit parameters. To tackle this, we rescale the standard devi- 
ations for the values of qqe such that we obtain a reduced 
chi-squared value of 1. 

To determine which variable to include first in our fit, 
we perform a 1-dimensional (ID) linear fit of our recon- 
structed values of «ce versus each of the binary parameters. 
Specifically, we fit the function of the form 

log 10 QcE = eo + eilogQ, (12) 

where eo and ei are constants. 

Then, to determine if a linear fit - as opposed to a 
constant fit - to the data is actually warranted (i.e. we reject 
the null hypothesis J^ , ), we perform the F-test (see Section 
[Bl) . We choose a significance level of a = 0.05. The first 
variable which we include in our fit is the one which gives 
the smallest probability, 5* (J 7 ), that we will obtain an F- 
statistic smaller than the obtained value, T . 

Once we have added our first variable to the fit, we de- 
termine the next binary parameter to include, by performing 
both a constant and a ID fit to the residuals as a function of 
the remaining binary parameters. The next variable which 
we add is the one which gives the smallest value of 8P(J-). 
In principle, we can repeat this process until we have a fit 
in terms of all binary parameters. 

However, special care should be taken when performing 
regression-analysis; strictly speaking, the predictor variables 
should be independent of one another. We check whether 
there is any correlation between the binary parameters via a 
correlation matrix plot as shown in Fig. [2] There are clearly 
strong correlations between Mwd, Mi and PcE.i- This is to 
be expected because the core mass of a Roche lobe-filling 
star is a (weak) function of the total mass, and as a func- 
tion of the orbital period via the star's Roche lobe radius 
and Kepler's Third Law. Hence, by including some or all of 
these variables we may over-fit the data, and we may render 
certain variables redundant in the fit. 

As we add each variable to our fit as described above, 
we once again apply the .F-Test to determine whether the fit 
is a statistically better fit than the one without this added 
binary parameter. 

Our results for the case with A = A g and A = Ab are 
discussed in turn in the following Sections. 



3 RESULTS 

The observed binary parameters (component masses and or- 
bital period) as well as the cooling age of the white dwarf, 
tcool) are summarised in table [Cll for PCEBs detected by 
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the S PSS. Table [C2l is the same but for PCEBs recorded 
in the lRitter fc Kolbl (|2003h catalogue, Edition 7.14 (2010). 
The average values of qce, Pce,; and Mi and their standard 
deviations (calculated from distributions similar to those il- 
lustrated in Fig. [TJ for each observed PCEB are also sum- 
marised in Table IC1I and IC2I where we do not consider the 
internal energy of the progenitor primary envelope. Tables 
IC3I and IC4I give the average values of Mi, PcE.i and «ce 
for each observed PCEB when we do consider the internal 
energy of the progenitor primary envelope. 



can reject for both of these alternative fits. Hence, we 
suggest 

Mw 



log 10 a C E = -0.31 - (0.95 ± 0.13)log 10 

+ (0.28±0.02)log 10 (^p) 
-(0.17±0.09)log 10 (j^. (15) 
Fig. [5] is the same as Fig. [4] except equation (|15[) is used. 



3.1 For A = A g 

Our values of ace versus Mwd, Mi, Mi, Pce and PcE,i 
are shown in Fig. [3] for the case where we consider A = A g . 
Each point corresponds to the mean value of oce, while the 
vertical lines correspond to the standard deviation. 

The values of eo, £i for each panel of Fig. [3] are given in 
Table Q] The solid green curves in Fig. [3] are the linear fits. 
For comparison, we also fit a constant function to the data 
in Figs. [3] of the form 



1o Sio q ce = eo, 



(13) 



where eo is a constant. For the data in Fig. we find that 
e = -0.06 ± 0.04. 

Values of T [and the probability that we exceed this 
value, £?(T)\ for the data in Fig. [3] are given in Table [T] Ta- 
ble[T]also indicates whether we can reject the null hypothesis 
Jfif (see Section [Bj. 



3.1.1 A fit in terms of the progenitor parameters 

Table Q] indicates that we can reject for PcE,i and M 2 . 
The first variable which we add to our fit is PcE.i as this has 
the smallest value of ^(T). Subsequent variables which we 
add, and the corresponding fit parameters, are summarised 
in the top panel of Table [5] Specifically, Table [2] gives the 
values f3o... m , the F-statistic J- and 8?{J-). We also indicate 
whether we can reject J% for each fit. We adopt a signifi- 
cance level of 0.05. 

We find that adding binary parameters in addition to 
fcE,i do not provide an improved fit over one in terms of 
PcE.i only. Indeed, S^i^F) > 0.05 in all cases. For such fits, 
we therefore cannot reject Jtf^. Hence, we find the following 
fit of acE in terms of PcE.i which is given by 

log 10 Q CE = (0.94 ± 0.33) - (0.34 ± 0.11)log 10 (^p) ■ (14) 

Our reconstructed values of qce versus values calculated 
from equation (I14[) are shown in Fig. U 



3.1.2 A fit in terms of the PCEB parameters 

The first variable which we add to our fit is Pce- Subsequent 
variables that we add, and the corresponding fit parameters, 
are summarised in the bottom panel of Table [2] In contrast 
to our fit of qce in terms of the progenitor parameters, we 
find that adding A/wd and then M2 is an improvement over 
the fit which is only in terms of Pce; we have T = 16.6 and 
38.4 for the former and latter cases respectively. Hence, we 



3.2 For A = A b 

The reconstructed values of «ce versus each binary param- 
eter are shown in the corresponding panels in Fig. [6] The 
values of eo, ei are summarised in Table [3] as are the T- 
statistics, and whether we can reject . For comparison 
with our linear fits, we once again perform a constant fit 
to our values of oce in the form of eqn. (|13[) .We find that 
e = -0.49 ± 0.05. 



3.2.1 A fit in terms of the progenitor parameters 

Table|3]indicates that we can reject J^q 7 " for the ID fit of qce 
versus Mwd, Mi, M2 or PcE,i- We obtain the smallest value 
of 2?(T) for Mwd, and so this is the first binary parameter 
we add to our fit. Subsequent additions and the correspond- 
ing fit parameters are summarised in the top panel of Table 
[4] which shows the same quantities as displayed in Table [2] 
Adding more binary parameters to our fit does not pro- 
vide an improvement over the fit only in terms of Mwd. 
This is probably a result of the fact that trends exist be- 
tween Mi, Mwd and PcE,i- Therefore, we suggest a fit of 
qce as a function of A/wd only, i.e. 



log 10 Q C E = (-1.03±0.10)-(2.10±0.34)log 10 



M v 



M . 



.(16) 



Fig. [7] is the same as Fig. [4] except equation (|16[) is used. 



3.2.2 A fit in terms of the PCEB parameters 

We once again start with Mwd in our fit, with the fit pa- 
rameters for additional variables summarised in the bottom 
panel of Table 2] Including Pce and then M2 provides an 
improved fit over one in terms of A/wd only; we can reject 
J^q 7 " in either case. As for the case where A = A g , we suggest 
a fit of qce as a function of A/wd, M2 and Pce, given by 



>'W< i- = - J -A*- 12.50 zztUO.I..^,, ( ^ 



WD 





+(0.31±0.02)log 10 (^p) 



-(1.02±0.09)log 10 



Ma 

Mr, 



(17) 



Fig. [8] is the same as Fig. [4] except equation (|17[) is used. 
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Table 3. Same as Table [T] except for the data shown in Fig. [6] 



Q 


£0 




T 


<?{T) 


Reject ? a 


%D 


-1.03±0.10 


-2.10±0.34 


36.1 


< 10" 6 


Yes 


Mi 


-0.21±0.07 


-1.26±0.23 


30.0 


< io- 6 


Yes 


M 2 


-0.89±0.13 


-0.70±0.20 


2.9 


6 x 10~ 4 


Yes 


PCE 


-0.45±0.06 


0.15±0.10 


2.7 


1 x lO" 1 


No 




1.55±0.43 


-0.76±0.16 


22.4 


2 x 10- 6 


Yes 



a -^crit = 4.00 for a 0.05 significance level. 



Table 4. The 


same as 


Table [2] except 


now for our data calculated using Aj,. 








Variables 


Po 


Pi 


02 Pz 04 


T 




Reject ? 








Progenitor parameters 








Af WD + M 2 


-1.08 


-1.94±0.29 


-0.14±0.14 


0.4 


0.5 


No 




-0.54 


-1.65±0.39 


-0.13±0.15 -0.17±0.16 


0.6 


0.4 


No 




0.98 


0.22±1.26 


-0.07±0.15 -0.46±0.24 -0.98±0.63 


1.1 


0.3 


No 








Observable (PCEB) parameters 








M W D + PCE 


-1.18 


-3.20±0.26 


0.22±0.03 


71.8 


< io~ 6 


Yes 


+M 2 


-1.53 


-2.50±0.16 


0.31±0.02 -1.02±0.09 


127.5 


< 10" 6 


Yes 



4 DISCUSSION 

We have found statistical evidence that the CE ejection effi- 
ciency has at least a first order dependence on the white 
dwarf mass, the progenitor primary mass, the secondary 
mass or the orbital period at the point the CE phase com- 
mences if we consider the case where A = Ab- For A = A g , 
we find correlations between qce and M2, Pce or PcE.i- 
However, whether we consider A = A g or A = Ab, we do not 
obtain a statistically significant fit of qqe in terms of more 
than one progenitor binary parameter. 

This behaviour may be a result of the large uncertain- 
ties of the progenitor system parameters for a given observed 
PCEB. As Fig.[T]demonstrates, an uncertainty for the white 
dwarf mass of a fewxO.01 Mq can translate into possible 
progenitor primary masses that span between approximately 
1 and 2 Mq , or into a pre-CE orbital period that can range 
from a few hundred to approximately 1000 d. As Tables [Cll 
to lC4l show. this can result in a standard deviation in qce 
as large as approximately 100 per cent of the mean value. 
This presents a challenge to efforts aimed at discerning po- 
tentially small, albeit real, trends with binary parameters. 

Even if the trends which we see are statistically sig- 
nificant, it is unclear whether they have a physical under- 
pinning. The initial distribution functions which we use to 
weight the possible values of qce for a given PCEB have 
been inferred from binary populations which may be prone 
to selection effects. Hence, for a given observed system, the 
calculated average value of cvce may not be the true average 

QCE- 

We now address whether there is a source of energy 
in addition to gravitational energy which is used during the 
ejection of the CE. We find that if we do consider the internal 
energy of the progenitor primary envelope, then this does 
bring the ejection efficiencies wi thin the canonical ran ge of 
< qce < 1, i n agreement withlZorotovic et all (|201Cf ) and 
as predicted bv lWebbinkl (|200&i ). However, we still find that 
«ce ^ 1 for PCEBs with brown dwarf secondaries (Ma < 0.1 

There is some debate as to what extent the internal 



energy of the giant 's envelope pla ys a role in the ejection 
of the CE. Indeed, lHarpaj (1 19981 ) argues that during the 
planetary nebula phase, the opacity of the giant's hydrogen 
ionization zone decreases during recombination. Hence, the 
envelope becomes transparent to its own radiation. The ra- 
diation will therefore freely escape rather t han push against 
the material to eject it. On the other hand, Wcb binkl (|2008h 
argues that the ionization zone of giant and AGB stars is 
buried beneath a region of the envelope whose opacity is 
dominated by heavy elements. Once recombination is trig- 
gered at the onset of the CE phase, the escaping radiation 
can therefore 'push' against these high-opacity layers. 

Ilvanova fc Chaichenetsl l|201lf ). on the other, suggest 
that the enthalpy of the primary giant's envelope, as op- 
posed to the internal energy of the gas, should be considered. 
Indeed, they found that there exists a region in the giant 
interior, coined the 'boiling pot zone', where the enthalpy 
is greater than the gravitational potential. This excess of 
energy may cause an outflow of material lying above the 
boiling-pot zone, and hence facilitate the ejection of the CE 
from the system. They find that by considering the enthalpy 
of the progenitor giant permits the formation of low- mass X- 
ray binaries with a low-mass companion, without resorting 
to qce > L 

Our calculations with A = Ab assume that all of the 
internal energy of the envelope is used to unbind the CE 
from the system, i.e a t h = 1 in eqn. ([5]), which is clearly 
unrealistic. Nonetheless, by reconstructing the CE phase 
with and without consideration for the internal energy of 
the primary envelope, we provide upper and lower limits 
for qce- Furthermore, our reconstructed values of acE and 
their trends with binary parameters presented in Section [3] 
will provide useful constraints on detailed models of the CE 
microphysics. 

Another result to note from our investigation regards 
the progenitors of those PCEBs with high mass white 
dwarfs. For PCEBs with white dwarf masses Mwd 0.7 
Mq, there are possible progenitor primaries which have 
Ab < at the onset of the CE phase (see also Section f2 . 3 [) . 
The envelope is therefore unbound from the giant as a result 
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Figure 6. Similar to Fig. [3] except for the case where we take A = Ab. The solid blue lines indicate linear fits to the data using the 
bootstrapping technique (see Section [4] for further details). 



of radiation pressure driven winds. A particularly striking 
example is IK Peg, where a vast majority of the possible 
progenitors have Ab < (see Fig. 0. This begs the ques- 
tion: d id IK Peg really form via a CE phase? iDermine et al.l 
(2009) investigated the effect that radiation pressure from 
giant stars has on the Roche lobe geometry. If / is the ra- 
tio of the radiation to the gravitational force, they found 
that the critical Roche lobe around the donor star becomes 
meaningless if / = 1. Hence, mass lost from the star will 
be in the wind regime, and material may not necessarily 
be channelled directly into the Roche lobe of the secondary 
star. Hence, in the context of giant stars, a situation where 
/ > 1 may mean that a CE phase is avoided. Theoretical 
studies investigating this possibility need to be carried out 
further. 

In the following Sections, we discuss important points 
regarding the reconstructed values of qce versus each binary 
parameter in turn. 



4.1 acE versus Mwd 

The following argument may further help to explain why 
we can reject more frequently for the data calculated 
using A = Ab rather than when A = A g is used. Recall that 
there are some PCEB progenitor primaries with Mwd > 0.7 
Mq which have Ab < (see Section Q. Furthermore, the 



progenitor primaries of these PCEBs which do have Ab > 
have very large values of Ab (and a correspondingly small 
binding energies) owing to their large radii when they fill 
their Roche lobes at the start of the CE phase. As a result, 
qce will be very small (for example IK Peg, where we find 
a mean CE ejection efficiency of qce ~ 0.19). This explains 
why the values of |ei| are larger when we consider A = Ab 
than when we consider A = A g . 

Furthermore, we find that as the white dwarf mass in- 
creases, the number of possible progenitor primaries with 
Ab > decreases. This is graphically illustrated in Fig. [5] 
which shows the two-dimensional distribution in Mi and 
P C E,i for three systems: SDSS2214-0103, SDSS0303-0054 
and IK Peg. Black points indicate that Ab > for the asso- 
ciated progenitor, while red points indicate that Ab < for 
that progenitor. 

This means that the range of possible progenitor pri- 
mary masses, and orbital periods at the start of the CE 
phase, and consequently the range in acE will decrease. As 
a result, the standard deviation in qce will decrease. Such 
data points will therefore have a large weighting in the fits. 
The linear fits between qce and Mwd may be greatly influ- 
enced by such data points. 

To determine if this is the case, and to determine if a 
linear relationship does exist among the other data points, 
we perform a bootstrap analysis to the qce versus Mwd 
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Intercept: - 1 .77401 02e-7 
Gradient: 0.99999979 




« CE (mod.) 



Figure 7. Same as Fig. [4] except acE( m °d.) has been calculated 
from eqn. Jl6j. 




a CE (mod.) 



Figure 8. Same as Fig. [4] except acE(inod.) has been calculated 
from eqn. dlTI . 



data in Fig. [5] The mean value of ei and the standard error 
from this analysis is given in Table [5] The linear fit calcu- 
lated from our bootstrapping analysis are shown in Fig. [6] 
as solid blue curves. 

In comparison to the linear fit, we find that the boot- 
strap analysis also indicates that cvce decreases with increas- 
ing Mwd, albeit a slightly steeper (i.e. a more negative) 
value of ei is predicted. In other words, a marginally stronger 
linear relationship between acE and Mwd is predicted by 
the bootstrapping analysis. By comparing the green and 
blue lines in Fig. [6] we can see that the data points at large 
white dwarf masses have the effect of flattening the gradient. 
The majority of the data points, however, suggest a steeper 
gradient. 



&. 3 
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Figure 9. Two dimensional distributions in primary mass, Mi 
and orbital period, Pce ii a * the very start of the CE phase, for 
three systems as indicated in the bottom right of each panel. Black 
points indicate that Ab > for those progenitors, while red data 
points indicate that Ab < 0. 



Table 5. The mean values with errors of the coefficient ei calcu- 
lated from our bootstrapping analysis of the data in each panel 
of Fig. E] 



Q 61 


Mwd 


-2.88±0.33 


Mi 


-1.87±0.28 


Ah 


-0.87±0.24 


Pce 


0.34±0.11 


PcE,i 


-0.79±0.15 



4.2 «ce versus Mi 

The linear fits to the corresponding panels in Fig. [3] and 
[6] suggest tha t oce decreases with increasing Mi . This is 
in contrast to lDe Marco et all (|201 lh . who found that oce 



increases with Mi. We explore this further in Section [5.21 

In light of the influence that data points corresponding 
to PCEBs with high mass white dwarfs may have on the 
fitting (see Section I4.1J1 we also perform a bootstrapping 
analysis on the Mi versus acE data in Fig. [6] (We repeat 
this analysis for the subsequent binary parameters in our 
discussion). Our bootstrap analysis predicts a slightly more 
negative value than the value obtained from the chi-squared 
fit. Therefore, those data points corresponding to high mass 
white dwarfs do not significantly affect the linear fit. 



4.3 acE versus M2 

We find that «ce decreases with increasing A'h , irrespective 
of whether we use A g or Ab. This finding is in agreement with 



Ii5e~ Marco et al.l (l20l"ll), bu t in co ntrast to the suggestion 
made by Politano fc Weilerl |2007) , who proposed that oce 
may increase with increasing M2. 
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This suggestion was motivated by the fact that very 
few, if any, CVs with brown dwarf secondaries at orbital pe- 
riods below 77 mi n have been de tected. This appeared to 
be in conflict with IPolitanol |2004h who estimated from his 
models that such systems should make up a pproximately 15 
per c e nt of the total C V population (see also lKolb fe Baraffd 
1 19991 ) . IPolitanol i|2004l ) suggested that this discrepancy may 
be a result of the decreasing energy dissipation rate of or- 
bital energy within the CE for decreasing secondary mass, 
and that below some cut-off mass, a CE merger would be 
unavoidable. Indeed, the small surface area of the secondary 
star will generate a correspondingly small friction force with 
the CE material. Also, as can be seen from eqn. (|6}, the 
change in orbital energy during spiral-in will be small for 
small secondary masses. 

However, as PCEBs with M2<0.1 M© are observed, 
the CE ejection must be very efficient i f such systems are to 
avoid a merger. |Pe Marco et al.l l|201ll ) suggested this may 
be achieved if smaller mass components take longer than the 
giant's dynamical time to penetrate into the envelope. The 
giant star would therefore have time to thermally re-adjust 
its structure, and use this thermal energy to aid the ejection 
of the CE. Note, however, that systems with secondaries 
M2 ss 0.1 Mq (see Fig. [5]), require a mean ejection efficiency 
close to 1, even if we consider the internal energy of the 
progenitor primary's envelope. 

The bootstrap analysis yields only a marginally steeper 
gradient than that obtained from the chi-squared fit. 



4.4 «ce versus Pqe 

We find that oce increases with increasing Pce , irrespective 
of whether we consider the internal energy of the progenitor 
primaries or not. 

Furthermore, the orbital separation at the point the CE 
terminates, Aqej (found via Pce), for the observed PCEBs 
indicates how deep the secondary star penetrated into the 
envelope of the giant star by the time the CE ceased. This 
may shed further light into the role that the structure of 
the gi ant primary plays during the termination of the CE 
phase. lYorke et al.l 1 19951) suggested that a steep decrease in 
density within the giant star is required for the successful 
termination of the CE phase. Such a decrease is character- 
ized by a flat mass-radius profile of the star. 

Once the secondary star reaches this region, the energy 
dissipation rate within the CE decreases, and the CE phase 
comes to an end, provided that the envelope material ex- 
terior to this region h as been unbound from the system. 
iTerman fc Taaml l|l996l ) suggest that the values of AcE,f for 
PCEBs may be estimated from the location within the CE 
where the mass-radius profile is no longer flat, i.e. when the 
quantity V = dlnP/dlnr is at a minimum. Here, P is the 
pres sure and r is the r adial coordinate. 

ITerman fe T aam ( 199(|) found from their hydrodynam- 
ical simulations that the orbital separation of PCEBs im- 
mediately after the CE phase was proportional to the white 
dwarf mass (see their Fig. 8), as a consequence of the fact 
that the point where V is minimum moves further away from 
the stellar centre for more evolved stars. Therefore, a plot 
of Mwd versus Ace,{ for our observed sample of PCEBs 
would provide an excellent comparison with the aforemen- 
tioned theoretical study. 




Figure 10. A plot of white dwarf mass versus the post-CE or- 
bital separation for our observed sample of PCEBs. The green 
dotted line is a constant fit to the data, while the green solid line 
is a linear fit to the data. The thick red border outlines the ap- 
proximate region in Aq e.I ~ tfwp space inhab ited by the data 
points given in Fig. 8 of ITerman fc Taaml <|l996h . calculated from 
their hydrodynamical simulations. 



Such a plot is shown in Fig. [10] To determine if a linear 
relationship between Mwd and Ace,! does exist, we fit a 
linear function given by 



log 10 A C E,t = eo + eilog 10 Af W D, 



(18) 



where Ace,{ is expressed in solar radii, and Mwd is ex- 
pressed in solar masses. We find eo = 1.07 ± 0.01 and 
ei = 2.76 ± 0.03. For our value of the number of degrees 
of freedom, vu nl = 58, we therefore obtain a reduced chi- 
squared value of x?mr = 98.62. This curve is shown as 
the green solid line in Fig. 1101 For comparison, we also 
show a constant fit to the data, log 10 AcE,f = eo, where 
e = 0.435 ± 0.004, and Xconst = 213.39 (reduced), for 
^const = 59 degrees of freedom. 

The chi-squared values indicate that a linear function 
provides a better fit to the data in Fig. [10] than a constant 
one. To determine whether an extra parameter is warranted, 
we once again apply the F-Test. For our values of vu nl and 
i/ colls , we obtain T = 67.5 with £?(T) < 10 -6 . Hence, be- 
cause £?(T) < a, we can reject . We therefore find 
intriguing statistical evidence that there is a functional de- 
pendence between Mwd and Acej- 

The thick red border in Fig. [TO] indicates the region 
in the Ace. i — Mwd plane i nhabit ed by the data points 
in Fig. 8 of ITerman fe Taaml l|l996l) . calculated from their 
hydrodynamical simulations. It can be seen that there is 
reasonable ov erlap between the observ ed PCEB sample and 
the results of ITerman fc Taaml (|l996l ) for Mwd < 0.7 M . 
In contrast, their simulations over-estimate the final orbital 
separation by as large as a factor of about 3 for systems with 
more massive white dwarf primaries. 

However, as noted by lYorke et al.l l|l995l ). the orbital 
separation when the CE terminates does not necessarily co- 
incide exactly with the point when V becomes a minimum. 
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lYorke et al] (|l995l ) and iTerman fc Taaml (|l996l ) estimate 
that the orbital separation at which the orbital decay suffi- 
ciently decreases is approximately 3 to 10 t imes smaller than 
the di stance when V becomes a minimum. ITerman Sz Taaml 
(1996) assume a factor of one sixth of this distance in their 
Fig. 8, although it is more likely that this factor will differ 

from system to sys tem. 

More recently, llvanoval (|2011T ) discussed the idea of the 
'compression point' within the primary giant, and the role 
that this may play during the CE ejection phase. The com- 
pression point is the location within the giant primary star 
where the value P/p (her e P is the pre ssure and p is the 
density) is at a maximum. Ilvanov al <|201ll ) found that, if the 
companion star is sufficiently massive, the envelope will be 
unbound up to the mass coordinate corresponding to the 
location of the compression point. Hence, the final orbital 
separation at the end of the CE phase may correspond with 
the radial coordinate associated with the compression point. 
Hence, our fit here may provide a useful diagnostic of the 
termination of the CE phase. 

4.5 «ce versus PcE.i 

The data points at high values of logio(PcE,i/d) in Fig. [6] 
may overly influence the linear fits. The reason is the same 
as for those data points at high white dwarf masses, as dis- 
cussed in Section 14.11 primaries that fill their Roche lobes 
towards longer orbital periods will have very large or neg- 
ative values of At,. Fewer progenitor systems will therefore 
have Ab > 0, decreasing the range in possible values of acE- 
We also find that the difference between ei calculated 
from the chi-squared fit and our bootstrap calculation is neg- 
ligible; our bootstrap analysis gives a slightly steeper value 
of €\ compared with that obtained from the chi-squared fit. 




0.1 1 10 

a CE (Zorotovic et al. 2010) 



Figure 11. A comparison of o ur reconstructed value s of oce 
(y-axis) with those calculated bv lZorotovic et alj l|2010h (x-axis). 
PCEBs formed only from a primary progenitor on the first giant 
branch are indicated by the red points, while blue points indicate 
that the progenitor primary was on t he AGB. The solid bla ck line 
is of the form acE(this workl=anR l|Zorotovic et a.1 .1 l201oh . Fits 
to the red and blue data points are shown as the red and blue 
lines respectively, while the green solid line is a linear fit to all the 
data points. Here, and in Figs. [T2l and 113 1 error-bars are omitted 
for clarity. 



T 



5 COMPARISONS WITH PREVIOUS STUDIES 

5.1 Comparisons with [Zorotovic et al.l (|20T0h 

Fig. [TJJ compares our values of qqe calcu lated with A = 
Ab (y- axis) with those values calculated bv lZorotovic et all 
(2010) ( x-axis) . The blac k line indicates y = x. 

(|201Gr ) determine 



For a given PCEB, IZorotovic et al 



mean values of «ce separately from progenitors that com- 
menced the CE phase o n either the first g iant b ranch (FGB) 
or on the AGB. Hence. IZorotovic et al.l (|2010h provide two 
mean values of acE for some PCEBs. This is in contrast 
to the present investigation, where we calculate one mean 
value of qce for each PCEB, from all possible progenitors. 

Therefore, to place our comparison in Fig. [TTJ on an 
equal footing, we only compare values of acE for systems 
that only have FGB progenitors (shown as red points in 
Fig. Ill |) and for systems with only AGB progenitors (blue 
points). The red and blue solid lines are linear fits to the 
corresponding data points, while the green curve is a fit to 
all the data points. 

For the case of those systems which have AGB pro- 
genitors, we typically obtain larger values of oce than 
IZorotovic et~aLl i|2010l) , by no more than a factor of 2. For 
PCEBs with FGB progenitors, we typically get values of 
qce which are smaller by a factor of no more than about 3. 



1 - 



E— 



0. i 



a (DeMarco et al. 2011) 



Figure 12. Similar to Fig. Illl e xcept here we compare our values 
of oce with those calculated bv lDe Marco et al.l (I201lh (red data 
points). 
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While in t h is pre sent work we always take ath = 1, 
IZorotovic etafl i|2010l) assume that the efficiency of using 
the giant envelope's internal energy to eject the CE is equal 
to the efficiency of using the orbital energy. Their values of 
A are then calculated from this assumption. Furthermore, 
IZorotovic etafl (]2010h calculate their values of A using ana- 
lytical fit s to the values of i?bmdi calculated from the stellar 
models of iPols et al.l (| 1998T ) . This is in contrast to the tabu- 
lated values of Ab for ath = 1 used in the present investiga- 
tion. 

5.2 Comparisons with lDe Marco et al.l (|201lh 

We now compare our values of oce, ca l culate d using A = A g , 
with those found by |Pe Marco et al.l (|201lh in Fig. [12] for 
PCEBs that are common to both studies. 

Fig- HHindicates that we obtain larger values of «ce by 
between a factor of 2 and 3. Reasons for th is disparity may 
be as follows. Firstly, foe Marco et al] (|201ll ) use an alterna- 
tive expression for the envelope binding energy, given by the 
left-hand side of their eqn. 6 which, for a given progenitor 
primary giant, will give a smaller binding energy t han that 
given by the left hand side of the IWebbinkl (1 19841 ) formal- 
ism (c.f. our eqn. [3]l . Hence, for a given P CEB progenitor, 
a sma ller value of «ce will be given if the foe Marco et al.l 
l|201ll ) formalism is used. Secondly, our values of A g are cal- 
culated fro m both the mass and r adius of the progenitor 
giant, while foe Marco et al.l (|201 lh calculate their values 
purely from the mass of the progenitor on the ZAMS. Hence, 
their method does not take into account any variations of A 
due to the location of the star (of a given mass) along the 
FGB or AGB, although they find that such variations are 
smalL 

foe Marco et al] (|201 lh found that oce increases with 
increasing Mi but decreases with increasing Mi. Hence, 
they find that acE decreases with increasing mass ratio, 
q = M2/M1, from their PCEB sample. 

Fig. ll3l shows data points on the acE-Mi plane (top left 
panel), the acE-M2 plane (top right panel) and on the acE- 
q plane (bottom middle panel). Blue data points have been 
calculated in this present study , whil e red points have been 
calculated by iDe Marco et al.l (|201ll ) , again for the afore- 
mentioned overlapping sample of PCEBs. The lines of best 
fit are shown as the blue and red lines respectively. 

The fit parameters to the data in each panel is sum- 
marised in Table |6] as are the (unreduced) chi-squared val- 
ues, xfinr an d Xconsi along with T and £?(T). We also indi- 
cate whether we can reject J ffyf . The top-left panel of Fig. 
1131 shows that, in contrast to IPe Marco et al] (|201lT l . a lin- 
ear fit to our data suggests a decreasing qce with increasing 
Mi. However, the data point indicated by the arrow (which 
corresponds to V471 Tau) may be overly influencing the fit. 
Indeed, the standard deviation in «ce for V471 Tau is the 
smallest in the sub-sample. A linear fit to the PCEB sub- 
sample without V471 Tau is given by the dotted blue line in 
the top-left pa nel of Fig. | 13[ Th i s is no w consistent with the 
fit obtained bv IPe Marco et all (|201ll ). However, as Table [6] 
shows, we cannot reject . 

On the other hand, the trends obtained from the lin- 
ear fits to the blue data points in the to p-right and bot- 
tom p anels are in line with those found bv IPe Marco et ah! 
(|201lT ) . We see that we obtain a very good agreement with 



|Pe Marco et al] l|201l[ ) for the value of ei in the acE-q plane. 
Here, we can also reject for both the red and blue data 
point on the oce — q plane. 

However, if we now perform a linear fit to our entire 
PCEB sample on the oce — q plane (shown as the black 
crosses in the bottom panel of Fig. [13] and the black line 
respectively), then a trend between oce and q becomes less 
convincing. Indeed, we can no longer reject ". It is there- 
fore u ncertain whether the trend seen by IPe Marco et al.l 
|201ll) is physical, or whether it is due to the small PCEB 
sample which they used. 



6 CONCLUSIONS 

We reconstruct the common envelope (CE) phase for the 
current sample of post-common envelope binaries (PCEBs) 
with observationally determined component masses and or- 
bital periods. We perform multi-regression analysis on our 
reconstructed values of the CE ejection efficiency, oce, in 
order to search for correlations between qce and the binary 
parameters (both the progenitor and observed parameters) . 
This analysis is carried out with and without consideration 
for the internal energy of the progenitor giants' envelope. 

If the internal energy of the progenitor giants' envelope 
is taken into account, a statistically significant correlation 
(in terms of the progenitor parameters) is found between 
qce and the white dwarf mass, Mwd, only. We find that 
qce decreases with increasing Mwd- 

If we do not consider the envelopes' internal energy, we 
find the most statistically significant correlation exists be- 
tween oce and the orbital period at the start of the CE 
phase, PcE.i- Specifically, we find that qce decreases with 
increasing PcE,i- In terms of the PCEB parameters, whether 
we consider the internal energy or not, we find a significant 
correlation between oce and Mwd, the orbital period im- 
mediately after the CE phase, Pce and the secondary mass, 
M2. While our values of qce, and the aforementioned corre- 
lations, cannot make definite physical claims regarding the 
CE phase, they will be nonetheless useful to constrain the 
micro-physics modelled in detailed CE simulations. 

We re-iterate that uncertainties in the measured white 
dwarf masses result in uncertainties in the progenitor binary 
parameters, and in turn the values of qce. Future investi- 
gations into reconstructing the CE phase would therefore 
greatly benefit from more accurate determinations of PCEB 
white dwarf masses, for example from eclipsing systems. 

If the internal energy of the progenitor primary envelope 
is considered, this brings the values of acE more in line 
with the canonical range of < acE < 1- However, we still 
require acE > 1 to account for PCEBs with brown dwarf 
secondaries, i.e. with M2<0.1 Mq. 

Finally, a large majority of possible progenitor primaries 
of IK Pegasi will have envelopes that have a positive bind- 
ing energy. We suggest that IK Peg may have avoided the 
CE phase, and instead mass was lost from the system via 
radiation pressure driven winds. 
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0.1 

q = M 2 /M, 

Figure 13. Reconstructed values of oqe versus Mi (top left), M2 (top right ) and the mass ratio q = M2/M1 (bottom middle). Blue 
data points correspond to our values, while red points have been calculated bv lDe Marco et al.l fcoill) . The black crosses in the bottom 
panel are our entire PCEB sample. The lines of best-fit are shown as the blue, red and black lines respectively. The dotted blue line in 
the top panel is a best fit to our data but with V471 Tau (indicated by the arrow) omitted from the fit. 



Table 6. Linear fit parameters for the data shown in the top-left, top-right and bottom panels of Fig. 1131 



Study 


£0 


£1 




Xcons 


T 


3*{T) 


Reject J0q~ 


ckce versus Mi 


Dc Marco ct al. (2011) 


-0.41± 0.09 


0.73± 0.42 


16.28 


19.32 


3.0 


0.1 


No 


This work (overlapping systems) 


0.32±0.15 


-1.16±0.40 


34.07 


52.09 


8.5 


0.01 


Yes 


This work (no V471 Tau) 


0.01±0.17 


0.56±0.74 


23.34 


24.24 


0.6 


0.5 


No 


ckce versus M2 


Dc Marco ct al. (2011) 


-0.59±0.22 


-0.40±0.30 


17.43 


19.32 


1.7 


0.2 


No 


This work (overlapping systems) 


-0.37±0.12 


-0.63±0.19 


30.28 


52.09 


11.5 


0.004 


Yes 


c«CE versus q 


Dc Marco ct al. (2011) 


-0.97±0.22 


-0.84±0.27 


12.00 


19.32 


9.8 


0.006 


Yes 


This work (overlapping systems) 


-0.73±0.27 


-0.82±0.31 


36.48 


52.09 


6.8 


0.02 


Yes 


This work (all systems) 


-0.28±0.16 


-0.24±0.17 


136.41 


141.22 


2.0 


0.2 


No 
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Finally, the initial mass ratio distribution lunction, g(qi), 
can be calculated from 

q{qi) = fiq", (A3) 
where fj, is a normalisation factor, and v is a free parameter. 

APPENDIX B: THE F-TEST 

The F-test determines how much an additional term has 
improved the value of the reduced chi-squared value. If the 
chi-squared values from fits using m and m + 1 parameters 
are Xm an d Xm+i respectively, then the F-statistic, T, is 
determined using 

2 _ 2 

X 2 m+1 /(N-m-l)' 

where N is the number of data points. We then calculate 
the probability, 8P(T), of exceeding the calculated value of 
T. 

We can reject the null hypothesis that the additional 
term is warranted, , if the value of £P(J-) is less than 
some significance level, a, i.e. SP{J-) < a. 



APPENDIX C: OBSERVED AND DERIVED 
PARAMETERS FOR THE OBSERVED SAMPLE 
OF PCEBS 



APPENDIX A: INITIAL DISTRIBUTION 
FUNCTIONS 

We use three standard distribution functions to calculate 
the formation probability of the zero-age main sequence bi- 
nary system. The in itial mass function (IMF) is given by 
i|Kroupa et al.lll993h 

Mi,;/M < 0.1, 

. 0.29056M" 1 - 3 0.1 < Mi,i/M < 0.5, 

J( 1<%) » 0.15571M" 2 - 2 0.5 < A/i,i/M < 1.0, 1 ' 

0.15571M" 2 - 7 1.0<Mi,i/M Q , 

The initial orbita l separation distribution (IOSD) function , 
h(a,i), is given by (|lben fc Tutukovll 19841 ; iHurlev et alJliobl ) 

/ as/Re < 3 or a ; /R > 10 6 , 

h[ai) ~ \ 0.078636a i - 1 3 < ai /R < 10 6 . {A2) 
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Table CI. A table showing the orbital parameters of the observed sample of PCEBs from the SDSS. Also given are the orbital 
periods immediately after the CE phase, PcEi the average progenitor primary mass, Mi, and the orbital period, PcE,i, at the start 
of the CE phase and the average CE ejection efficiency, crcEi calculated with A g . 



System 


A/wd/M q 


m 2 /m q 




; 


orb/d 


P CE /d 


tcooi/Gyr 


Mi/M Q 


log(P C E,i/d) 


£>CE 




Ref. 


SDSS1435+3733 





505±0 


025 





218±0 


028 





126 





134 


0.275 


1 


26±0 


30 


2 


81±0 


14 





26±0 


18 


1 


SDSS0052-0053 


1 


220±0 


370 





320±0 


060 





114 





148 


0.421 


5 


18±1 


02 


3 


13±0 


12 





46±0 


13 


2 


SDSS2123+0024 





310±0 


100 





200±0 


080 





149 





149 


0.000 


1 


41±0 


30 


1 


39±0 


50 


4 


57±4 


57 


3 


SDSS1529+0020 





400±0 


040 





260±0 


040 





165 





171 


0.300 


1 


32±0 


25 


2 


21±0 


18 





87±0 


46 


2 


SDSS1411+1028 





520±0 


110 





380±0 


070 





167 





178 


0.009 


1 


42±0 


42 


2 


69±0 


22 





32±0 


22 


3 


SDSS1548+4057 





646±0 


032 





174±0 


027 





185 





191 


0.416 


2 


46±0 


26 


2 


84±0 


16 





95±0 


33 


1 


SDSS0303-0054 





912±0 


034 





253±0 


029 





134 





208 


>2.24 


4 


27±0 


18 


3 


08±0 


10 





65±0 


13 


1 


SDSS2216+0102 





400±0 


060 





200±0 


080 





210 





213 


0.297 


1 


32±0 


25 


2 


18±0 


26 


1 


29±0 


77 


3 


SDSS1348+1834 





590±0 


040 





319±0 


060 





249 





252 


0.184 


2 


03±0 


37 


2 


83±0 


20 





66±0 


30 


3 


SDSS0238-0005 





590±0 


220 





380±0 


070 





212 





239 


0.045 


1 


49±0 


50 


2 


58±0 


31 





49±0 


35 


3 


SDSS2240-0935 





410±0 


080 





250±0 


120 





261 





265 


0.443 


1 


32±0 


25 


2 


19±0 


31 


1 


20±0 


83 


3 


SDSS1724+5620 





420±0 


010 





360±0 


070 





333 





333 


0.000 


1 


30±0 


25 


2 


39±0 


09 





81±0 


11 


2 


SDSS2132+0031 





380±0 


040 





320±0 


010 





222 





224 


0.179 


1 


34±0 


25 


2 


08±0 


19 


1 


10±0 


56 


3 


SDSS0110+1326 





470±0 


020 





310±0 


050 





333 





333 


0.051 


1 


21±0 


23 


2 


67±0 


10 





48±0 


25 


1 


SDSS1212-0123 





470±0 


010 





280±0 


020 





333 





334 


0.191 


1 


19±0 


22 


2 


70±0 


09 





48±0 


25 


4 


SDSS1731+6233 





450±0 


080 





320±0 


010 





268 





271 


0.228 


1 


28 ±0.26 


2 


48±0 


29 





62±0 


47 


3 


SDSS1047+0523 





380±0 


200 





260±0 


040 





382 





384 


0.417 


1 


38±0 


31 


1 


74±0 


71 


1 


68±4 


68 


5 


SDSS1143+0009 





620±0 


070 





320±0 


010 





386 





387 


0.138 


2 


14±0 


11 


2 


84±0 


19 





87±0 


38 


3 


SDSS2114-0103 





710±0 


100 





380±0 


070 





411 





411 


0.018 


2 


65±0 


38 


2 


95±0 


16 





72±0 


23 


3 


SDSS2120-0058 





610±0 


060 





320±0 


010 





449 





450 


0.156 


2 


11±0 


10 


2 


84±0 


19 





98±0 


43 


3 


SDSS1429+5759 


1 


040±0 


170 





380±0 


060 





545 





564 


0.401 


1 


91±0 


65 


3 


13±0 


11 





91±0 


21 


3 


SDSS1524+5040 





710±0 


070 





380±0 


060 





590 





596 


0.109 


2 


74±0 


30 


2 


97±0 


13 





91±0 


25 


3 


SDSS2339-0020 





840±0 


360 





320±0 


060 





655 





657 


0.508 


1 


97±0 


97 


2 


86±0 


18 





87±0 


58 


2 


SDSS1558+2642 


1 


070±0 


260 





319 ±0.060 





662 





665 


0.609 


1 


77±0 


95 


3 


12±0 


11 


1 


17±0 


30 


3 


SDSS1718+6101 





520±0 


090 





320±0.010 





673 





673 


0.075 


1 


44±0 


42 


2 


73±0 


19 





85±0 


57 


3 


SDSS1414-0132 





730±0 


200 





260±0.040 





728 





729 


0.329 


2 


25±0 


68 


2 


84±0 


20 


1 


55±0 


68 


5 


SDSS0246+0041 





900±0 


150 





380±0.010 





728 





737 


0.309 


3 


98±0 


62 


3 


10±0 


10 


1 


00±0 


21 


2 


SDSS1705+2109 





520±0 


050 





250±0.120 





815 





815 


0.023 


1 


40±0 


40 


2 


81±0 


16 





95±0 


73 


3 


SDSS1506-0120 





430±0 


130 





320±0.010 


1 


051 


1 


051 


0.251 


1 


36±0 


30 


2 


21 ±0.47 


2 


82±2 


82 


3 


SDSS1519+3536 





560±0 


040 





200±0.080 


1 


567 


1 


567 


0.065 


1 


78±0 


35 


2 


74±0.18 


3 


24±1 


64 


3 


SDSS1646+4223 





550±0 


090 





250±0.120 


1 


595 


1 


595 


0.093 


1 


55±0 


48 


2 


78±0.17 


1 


86±1 


33 


3 


SDSS0924+0024 





520±0 


050 





320±0.010 


2 


404 


2 


404 


0.059 


1 


44±0 


40 


2 


81±0.16 


1 


78±1 


31 


3 


SDSS2318-0935 





490±0 


060 





380±0.070 


2 


534 


2 


534 


0.026 


1 


31±0 


32 


2 


73±0.18 


1 


66±1 


15 


3 


SDSS1434+5335 





490±0 


030 





320±0.010 


1 


357 


1 


357 


0.030 


1 


14±0 


22 


2 


82±0.13 


1 


78±1 


11 


3 



(i; 
(i 



Pvrzas et al] d2009h: f 2 lRebassa-Mansergas et all d200Sf l: (3)Nebot-G6mez-Moran et al. (2011, submitted); 
Nebot Gomez-Moran et al] d2009h ; (5 ~ISchreiber et al] feOOiji 



Table C2. The same as tabic [Cl] but for the observed sample of PCEBs from the iRitter 
(2010) 
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Kolbl l|2003t) catalogue, Edition 7.14 



System 


A/ W d/M 


M 2 /M 


I 


orb/d 




tcool/Gyr 


Mi/M 


log(P C E,i/d) 


C*CE 


Ref. 


0137-3457 


0.390+0.035 


0.053+0.006 





080 


0.084 





179 


1.33+0.25 


2.04+0.17 


2.63+1.33 


1 


HR Cam 


0.410+0.010 


0.096+0.004 





103 


0.104 





118 


1.32+0.25 


2.23+0.09 


1.41+0.72 


2 


RR Cae 


0.440+0.023 


0.180+0.010 





304 


0.309 


2 


037 


1.27+0.24 


2.47+0.13 


1.07+0.59 


3,4 


DE CVn 


0.540+0.040 


0.410+0.050 





364 


0.486 





895 


1.48+0.42 


2.83+0.16 


0.49+0.36 


4,5 


J2130+4710 


0.554+0.017 


0.555+0.023 





521 


0.527 





088 


1.79+0.37 


2.84+0.19 


0.62+0.33 


1.6 


EG UMa 


0.630+0.050 


0.360+0.040 





668 


0.689 





313 


2.28+0.35 


2.88+0.19 


1.17+0.49 


7,8 


1857+5144 


0.610+0.040 


0.410+0.030 





266 


0.266 





000 


2.16+0.36 


2.87+0.20 


0.55+0.24 


9,10 


BPM 71214 


0.770+0.060 


0.540 : 





202 


0.289 





120 


3.17+0.27 


3.07+0.10 


0.38+0.08 


11,12,13 


QS Vir 


0.780+0.040 


0.430+0.040 





151 


0.315 





370 


3.26+0.19 


3.07+0.10 


0.48+0.10 


4,14 


V471 Tau 


0.840+0.050 


0.930+0.070 





521 


0.521 





001 


3.70+0.25 


3.15+0.10 


0.35+0.07 


15,16 


LM Com 


0.350+0.030 


0.170+0.020 





259 


0.260 





032 


1.36+0.26 


1.82+0.16 


2.88+1.33 


17 


MS Peg 


0.490+0.040 


0.190+0.020 





174 


0.174 





027 


1.25+0.28 


2.73+0.15 


0.40+0.26 


17 


GK Vir 


0.510+0.040 


0.100 : 





344 


0.344 





002 


1.39+0.36 


2.72+0.15 


1.29+0.92 


4,18 


NN Ser 


0.535+0.012 


0.110+0.004 





130 


0.130 





001 


1.68+0.27 


2.66+0.17 


1.02+0.49 


19 


1042-6902 


0.560+0.050 


0.140+0.010 





337 


0.340 





080 


1.57+0.45 


2.75+0.17 


1.10+0.81 


21 


J2013+4002 


0.560+0.030 


0.230+0.010 





706 


0.710 





002 


1.84+0.35 


2.76+0.19 


1.74+0.88 


20,21 


FS Cet 


0.570+0.030 


0.390+0.020 


1 


230 


4.230 





001 


1.90+0.37 


2.82+0.20 


3.80+2.01 


20,22 


IN CMa 


0.580+0.030 


0.430+0.030 


1 


260 


1.260 





002 


1.97+0.37 


2.85+0.20 


1.51+0.73 


11,20 


BE UMa 


0.590+0.070 


0.250+0.080 


2 


291 


2.291 





000 


1.88+0.41 


2.77+0.19 


3.55+1.80 


23,24 


UZ Sex 


0.680+0.230 


0.220+0.050 





597 


0.597 





084 


1.72+0.69 


2.78+0.18 


1.12+0.72 


25,26 


J1016-0520AB 


0.610+0.060 


0.150+0.020 





789 


0.790 





002 


2.13+0.39 


2.78+0.19 


2.88+1.26 


27 


2009+6216 


0.620+0.020 


0.189+0.004 





741 


0.741 





020 


2.30+0.30 


2.82+0.19 


2.24+0.93 


28 


CC Cet 


0.400+0.110 


0.180+0.050 





287 


0.287 





000 


1.37+0.30 


1.95+0.44 


2.57+2.49 


25,29 


HZ 9 


0.510+0.100 


0.280+0.040 





564 


0.564 





086 


1.42+0.40 


2.66+0.22 


0.93+0.64 


30,31 


LTT 560 


0.520+0.120 


0.190+0.050 





148 


0.162 


1 


040 


1.44+0.42 


2.61+0.24 


0.60+0.36 


32 


IK Peg 


1.190+0.050 


1.170 : 


21.722 


21.722 





027 


6.10+0.24 


3.28+0.12 


3.98+1.19 


33,34 



200711; (6lMaxted et ahl d2004h ; ( 7"lBleach et al 
flO lShimanskv et al.l j2009Tl; (ll jKawka et al " 
(14 ]Q'Donoehue et al] d2003ft: d5lO'Brien et al 



j2002h : (8 lShimanskii fe Borisovl J2002J); P lAungweroiwit et al 
12002?) ; (12 lKawka fe Vermes! d2003ft ; (13 lKawka fe Vennesl 
_ jpOOlll; (16lKaminski et al] d2007l) ; (1 7 lShimanskv et al.l 

181Fulbrieht et all 11993T); (19lCatalan et al.l dl994h;(20lKawka et al. d2008ft ; (2llGood et al] j2005ri ; (2 2lVennes fe Thorsterisenl 



1994ft; (23 



Shimanskii ct al 



d2008l); 



27lVennes et al.l [|l999h; (28 , _ . . _ 

1987ft ; (32 lTappert et all [|2007r i; (33 lLandsman et al] dl993h ; f34 lVennes et al] lll998l) 



(5lvan den Bcssclaar ct al. 

l|2007 ); 
(2005); 



(24lRaguzova et al] ll2003ft; (25~lSaffer et al] dl993ft~(26feruch fe Diad dl999ft; 
VIorales-Rueda et al] ll2005ft; C29 lSomers et al.ll 1996h; QO lLanning fe Pescbl l|l98lft : (3l ]staufferl 
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Table C3. Similar to table ICll but now showing the values of Mi, Pce i an( i a CE f° r the SDSS PCEBs, where consider the internal 
energy of the primary progenitor envelope, i.e. we use A = At,. 



System A/i/Mq log(P C E.i/d) «ce 



SDSS1435+3733 


1 


.26±0. 


.30 


2 


.81±0. 


.14 


0. 


17±0. 


05 


SDSS0052-0053 


5 


.35±1. 


.24 


3 


.01±0. 


.09 


0. 


03±0. 


.02 


SDSS2123+0024 


1 


.41±0. 


30 


1 


.39±0. 


.50 


2. 


04±2. 


.04 


SDSS1529+0020 


1 


.32±0. 


25 


2 


.21±0. 


.18 





37±0. 


21 


SDSS1411+1028 


1 


.42±0. 


.42 


2 


.69±0. 


.22 


0. 


14±0. 


07 


SDSS1548+4057 


2 


.46±0. 


26 


2 


.84±0. 


16 


0. 


13±0. 


11 


SDSS0303-0054 


4 


.37±0. 


16 


2 


.96±0. 


03 


0. 


03±0. 


.02 


SDSS2216+0102 


1 


.32±0. 


25 


2 


.18±0. 


26 





56±0. 


.34 


SDSS1348+1834 


2 


.03±0. 


37 


2 


.83±0. 


20 





14±0. 


.10 


SDSS0238-0005 


1 


.48±0. 


.49 


2 


.58±0. 


.31 


0. 


18±0. 


.14 


SDSS2240-0935 


1 


.32±0. 


25 


2 


.19±0. 


31 


0. 


54±0. 


35 


SDSS1724+5620 


1 


.30±0. 


25 


2 


.39±0. 


.09 





34±0. 


18 


SDSS2132+0031 


1 


.34±0. 


25 


2 


.08±0. 


.19 


0. 


.48±0. 


.24 


SDSS0110+1326 


1 


,21±0. 


23 


2 


.67±0. 


.10 


0. 


25±0. 


09 


SDSS1212-0123 


1 


.19±0. 


22 


2 


.70±0. 


.09 





26±0. 


08 


SDSS1731+6233 


1 


.28±0. 


26 


2 


.48±0. 


.29 


0. 


30±0. 


.18 


SDSS1047+0523 


1 


.38±0. 


31 


1 


.74±0. 


71 


2. 


08±2. 


08 


SDSS1143+0009 


2 


,14±0. 


11 


2 


.84±0. 


.19 


0. 


.17±0. 


13 


SDSS2114-0103 


2 


.62±0. 


37 


2 


.93±0. 


15 


0. 


07±0. 


.07 


SDSS2120-0058 


2 


.11±0. 


.40 


2 


.84±0. 


.19 


0. 


20±0. 


15 


SDSS1429+5759 


4 


.84±0. 


61 


2 


.99±0. 


03 





05±0. 


03 


SDSS1524+5040 


2 


.73±0. 


30 


2 


.96±0. 


12 





07±0. 


07 


SDSS2339-0020 


1 


.80±0. 


.72 


2 


.84±0. 


17 


0. 


27±0. 


.20 


SDSS1558+2642 


1 


.50±0. 


.83 


2 


.98±0. 


03 


0. 


06±0. 


.04 


SDSS1718+6101 


1 


.44±0. 


.42 


2 


.73±0. 


.19 





37±0. 


18 


SDSS1414-0132 


2 


.15±0. 


.58 


2 


.82±0. 


.19 


0. 


30±0. 


.28 


SDSS0246+0041 


3 


.87±0. 


60 


2 


.99±0. 


05 


0. 


05±0. 


03 


SDSS1705+2109 


1 


.45±0. 


.40 


2 


.79±0. 


16 





51±0. 


.20 


SDSS1506-0120 


1 


.36±0. 


30 


2 


.21±0. 


.47 


1. 


26±1. 


.24 


SDSS1519+3536 


1 


.78±0. 


35 


2 


.74±0. 


.18 


1. 


02±0. 


.54 


SDSS1646+4223 


1 


.55±0. 


,18 


2 


.78±0. 


17 


0. 


76±0. 


.38 


SDSS0924+0024 


1 


.44±0. 


.40 


2 


.81±0. 


16 


0. 


85±0. 


33 


SDSS2318-0935 


1 


.31±0. 


32 


2 


.73±0. 


18 





85±0. 


35 


SDSS1434+5335 


1 


.14±0. 


.22 


2 


.82±0. 


.13 


1. 


20±0. 


.39 
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Table C4. Same as tabic [C3] but now showing the values of Mi, PcE.i and «ce for the PCEBs from the lRitter fc Kolbl l|2003f ) catalogue, 
Edition 7.14 (2010). 



System Aii/M log(fcE,i/d) «ce 



0137-3457 


1 


.33+0 


.25 


2 


.04+0. 


.17 


1 


.12+0.57 


HR Cam 


1 


.32+0 


.25 


2 


.23+0 


.09 





.59+0.31 


RR Cae 


1 


.27+0 


.24 


2 


47+0. 


.12 





.48+0.23 


DE CVn 


1 


.48+0 


.42 


2 


.83+0 


.16 





.22+0.09 


J2130+4710 


1 


.79+0. 


.37 


2 


.84+0 


.19 





.18+0.09 


EG UMa 


2 


.28+0 


.35 


2 


.88+0 


.19 





.19+0.16 


1857+5144 


2 


.16+0 


.36 


2 


.87+0 


.20 





.10+0.08 


BPM 71214 


3 


.17+0. 


.26 


3 


.02+0 


.08 





.018+0.017 


QS Vir 


3 


.32+0 


.17 


3 


.00+0 


.06 





.022+0.018 


V471 Tau 


3 


.77+0. 


.24 


3 


.06+0 


.04 





.018+0.011 


LM Com 


1 


.36+0 


.26 


1 


.82+0. 


.16 


1 


.26+0.58 


MS Peg 


1 


.25+0. 


.28 


2 


.73+0. 


.15 





.23+0.07 


GK Vir 


1 


.39+0 


.36 


2 


.72+0. 


.15 





.68+0.23 


NN Ser 


1 


.68+0 


.27 


2 


.66+0 


.17 





.37+0.14 


1042-6902 


1 


.57+0 


.45 


2 


.75+0. 


.17 





.46+0.21 


J2013+4002 


1 


.84+0 


.35 


2 


.76+0 


.19 





.50+0.28 


FS Cet 


1 


.90+0 


.37 


2 


.82+0 


.20 


1 


.00+0.64 


IN CMa 


1 


.97+0. 


.37 


2 


.85+0 


.20 





.35+0.25 


J1016-0520AB 


2 


.13+0. 


.39 


2 


.78+0 


.19 





.60+0.44 


2009+6216 


2 


.30+0 


.30 


2 


.82+0 


.19 





.38+0.32 


CC Cet 


1 


.37+0. 


.30 


1 


.95+0 


.44 


1 


.15+1.11 


HZ 9 


1 


42+0. 


.10 


2 


.66+0 


.22 





.41+0.21 


LTT 560 


1 


44+0. 


.12 


2 


.61+0 


.24 





.25+0.14 


BE UMa 


1 


.88+0 


.41 


2 


.77+0. 


.19 





.97+0.63 


UZ Sex 


1 


.66+0 


.63 


2 


.77+0 


.17 





.39+0.27 


IK Peg 


6 


.29+0 


.17 


3 


.10+0 


.01 





.19+0.06 



